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ABSTRACT 

RNA thermometers are translational control 
elements that regulate the expression of bacterial 
heat shock and virulence genes. They fold into 
complex secondary structures that block translation 
at low temperatures. A temperature increase 
releases the ribosome binding site and thus 
permits translation initiation. In fourU-type RNA 
thermometers, the AGGA sequence of the SD 
region is paired with four consecutive uridines. We 
investigated the melting points of the wild-type and 
mutant sequences. It was decreased by 5 C when a 
stabilizing GC basepair was exchanged by an AU 
pair or increased by 11 C when an internal AG 
mismatch was converted to a GC pair, respectively. 
Stabilized or destabilized RNA structures are 
directly correlated with decreased or increased 
in vivo gene expression, respectively. Mg 2+ also 
affected the melting point of the fourU thermometer. 
Variations of the Mg 2+ concentration in the physio- 
logical range between 1 and 2 mM translated into a 
2.8 C shift of the melting point. Thus, Mg 2+ binding 
to the hairpin RNA is regulatory relevant. Applying 
three different NMR techniques, two Mg 2+ binding 
sites were found in the hairpin structure. One of 
these binding sites could be identified as outer 
sphere binding site that is located within the fourU 
motif. Binding of the two Mg 21 ions exhibits a 
positive cooperativity with a Hill coefficient of 1.47. 
Free energy values AG for Mg 2+ binding determined 
by NMR are in agreement with data determined from 
CD measurements. 

INTRODUCTION 

The maintenance of the intracellular metal ion concentra- 
tion is an indispensable requirement for the viability of 
eukaryotic and prokaryotic cells. Cells exhibit high 



intracellular K + and Mg 2+ concentrations while cytosolic 
Na + and especially Ca 2 concentrations are kept low (1). 
By varying the K + concentration, cells adopt the intracel- 
lular to the extracellular osmolality (2). K + acts as 
counterion for nucleic acids (DNA, RNA) (3,4) and has 
been found to stabilize the structure of cytosolic proteins. 
Mg 2+ is prominent for binding and activating ATP (5) as 
well as intracellular enzymes (1). In addition to K + , it is 
the most important cation acting as counterion of nucleic 
acids (1,4). It is involved in the maintenance of the struc- 
tural integrity of the plasma membrane (6) and stabilizes 
proteins. In gram-negative bacteria, Mg 2+ binding to lipo- 
polysaccharides stabilizes the periplasmic membrane (7). 
Most importantly, it is strictly required for rRNA folding, 
ribosome assembly and ribosome function (8,9). 

The importance of Mg 2+ ions for the viability of pro- 
karyotic cells has some important implications. First, 
Mg + starvation leads to growth arrest of bacterial 
cultures (10). Second, although the amount of Mg 2+ 
within a prokaryotic cell may vary by a factor of three, 
the RNA[kg]:Mg 2+ [kg] ratio («68) is almost constant 
under differing environmental conditions (10). 
Escherichia coli cells that enter the stationary phase 
reduce the concentration of ribosomes and thus reduce 
their Mg 2+ content (11). Third, E. coli cells react on 
Mg + starvation by the synthesis of polyamines including 
spermidine and putrescine which are able to substitute for 
Mg 2+ in many cases. For instance, the spermidine/Mg 2+ 
ratio bound to the ribosomal fraction of E. coli cells 
reduces with increasing extracellular Mg 2+ concentration 

(12) . Due to the diversity of different Mg 2+ binding part- 
ners and the multitude of different intracellular Mg 2+ 
binding sites, most of the intracellular Mg 2+ is present in 
a bound form, albeit transiently. E. coli cells in the expo- 
nential growth phase have a Mg 2+ content of 90-1 10 mM 

(13) , while only a small fraction of 1-2 mM is present in 
the physiologically active free form (14). In addition, most 
of the bound Mg 2+ fraction is in fast exchange with the 
free Mg 2+ in the cytosol (14,15) and a 20-30% loss of 
Mg 2+ is sufficient for cell-growth arrest (14). Due to this 
rapid exchange between free and bound Mg 2+ , also the 
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concentration of the free form has to be controlled tightly 
by the cell itself to assure cell viability. Indeed, Mg 2+ ions 
have been found to be actively transported among others 
in E. coli (15-19), Bacillus subtilis (20) and Salmonella 
typhimurium (21) and the Mg 2+ transport mechanisms of 
S. typhimurium and E. coli are equivalent to each other 
(1,14,18,21-25). 

Ions are important for stabilizing RNA tertiary struc- 
ture (3,9,26). Most often, K + and Mg~ + cations serve this 
task although Mn~ + ions and polyamines such as spermi- 
dine have the potential to substitute for some of the K + 
and Mg 2+ binding sites. Mn 2+ ions are able to substitute 
for 76% of the bound Mg 2+ ions in E. coli without im- 
pairing viability of the cells (18). These findings provide 
support for using Mn 2+ ions to map Mg 2+ binding sites 
within an RNA molecule (27-30). 

Due to its small size and divalent charge, the Mg 2+ ion 
has a notably high charge density. Therefore, the first 
water hydration sphere is more stable than for any other 
physiologically relevant cation (31). Mg 2+ coordinates six 
water molecules with octahedral geometry. Mg 2+ ions are 
able to bind to RNA with different binding modes. 
Generally, 'inner sphere' and 'outer sphere' contacts can 
be distinguished (9). In an 'outer sphere' complex, water 
molecules of the Mg(H 2 0) 2+ 6 complex but not the Mg 2+ 
ion itself make contacts to the functional groups of the 
RNA. In contrast, inner sphere contacts are characterized 
by at least one direct contact of the Mg 2+ ion to functional 
groups of the RNA. Inner sphere contacts can be classified 
by the number of inner sphere ligands other than H 2 0 
(type I, type II, type III, etc.), while 'type 0' is identical 
to an outer sphere contact (9). Klein et al. (9) modeled 
Mg 2+ ions into the structure of the large ribosomal 
subunit from Haloarcula marismortui and found that 
65% of the bound Mg 2+ ions are outer sphere 
coordinated. Within the remaining inner sphere contacts, 
most contacts belong to the classes type I or type Ha 
where inner sphere ligands are orthogonally oriented to 
each other. Most of these Mg~ + ions (type I and type 
Ha) are coordinated to either one or both non-bridging 
oxygens of the phosphate backbone. Outer sphere 
coordinated Mg 2+ ions are most often located within the 
major groove of the RNA helices and make indirect, 
water-mediated contacts to the heteroatoms of the 
nucleobases (9,27,32). 

Mg 2+ has a profound effect on the structure and 
function of riboswitches (33-35). Riboswitches are biosen- 
sors that control gene expression in response to cellular 
metabolites that they bind with high affinity and specifi- 
city (36,37). RNA thermometers are capable of sensing 
temperature. They control translation initiation of heat 
shock and virulence genes by sequestration of the Shine- 
Dalgarno (SD) sequence at low temperatures (38). At 
higher temperatures, the hairpin structure melts in a 
zipper-like manner allowing translation initiation to 
occur (39). Helix stability in simple synthetic RNA therm- 
ometers is directly correlated to gene expression levels of 
the downstream gene (40). 

We recently investigated the structure, function and 
melting mechanism of the Salmonella fourU RNA 
thermometer, which controls expression of the heat 



shock gene agsA (39,41). In the present work, we 
examined the in vitro and in vivo effects of a point 
mutation that weakens an extraordinarily stable GC 
base pair, which had been shown to prevent unfolding of 
the wild-type RNA at lower temperatures (39) (Figure 2). 
Furthermore, we highlight the Mg 2+ binding properties of 
the fourU RNA thermometer. Using NMR spectroscopy, 
we identified and localized two Mg 2+ binding sites and 
determined their dissociation constants as well as 
the cooperativity between these binding sites. Mg 2+ - 
dependent changes in thermal stability of the fourU 
RNA thermometer were determined by CD spectroscopy. 
Our study provides first evidence for Mg 2+ binding of an 
RNA thermometer and its relevance for RNA thermom- 
eter function. 



MATERIALS AND METHODS 

NMR spectroscopy 

NMR measurements were performed on a 600 MHz 
(14.09 T) Bruker NMR spectrometer equipped with a 
5-mm xyz-gradient probe and a 900-MHz (21.14T) 
Bruker NMR spectrometer equipped with a 5-mm cryo- 
genic HCN z-gradient probe. Spectra were recorded and 
processed using the Bruker TopSpin 2.0/2.1 software. 
Processed spectra were analyzed with Sparky 3.1.1.3. 

MgCb titration experiments 

Mg 2+ titration experiments of the 4U-hp2-wt RNA were 
performed in NMR buffer (see RNA preparation section) 
at a temperature of 10°C. ['H- 15 N]-HSQC experiments 
were recorded at a B 0 field strength of 900 MHz (21.14 
T) using a standard pulse sequence (42) (Topspin 2.1 
pulse program library: fhsqcf3gpph) with a hard 
Watergate pulse train for water suppression (43). Spectra 
were recorded with a spectral width of 30ppm in the 
15 N-dimension and 22ppm in the 'H-dimension. A total 
of 240 increments were recorded in the 15 N-dimension and 
4096 points in the 'H-dimension. ['H,' 5 N]-HSQC spectra 
were recorded at the following MgCl 2 concentrations: 0, 
0.3, 0.5, 0.7, 1.0, 1.5, 2.0, 3, 5 and 7mM. 

Overall chemical shift perturbations (CSPs) were 
calculated from the CSPs of the 15 N and the 
'H-dimension (44) according to the following equation: 



CSP 



/(CSP„) 2 +(CSP w /5) 2 



(1) 



Since association and dissociation of the Mg 2+ ions from 
the target RNA are fast on the NMR timescale, CSPs can 
be related to the Mg 2+ concentration within the sample 
according to the following equation: 



Ax) 



(x+l+a) - J((x+l+af-4x) 



(2) 



Equation (2) is derived from the mass action law for an 
RNA/Mg 2+ equilibrium with the RNA having one Mg 2+ 
binding site. In Equation (2)f(x) is the CSP of an individ- 
ual imino signal at a given x = [Mg 2+ ]/[RNA] ratio, a is 
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macroscopic 




K2 - K 2 ,4 + K 3i4 



Ki = 



Ki,2 * Ki i3 



Ki,2 + K13 



microscopic 




K* IX —W * If 
2,4 *1,2 - 1^3,4 ^1,3 



Figure 1. Mg binding model assuming two distinct Mg binding sites on the RNA. The RNA can adopt four different states (I, II, III, IV). 
Transitions between these states are described by microscopic dissociation constants (K 12 , ifi,3, K 2 4, ^3,4). Transitions between free RNA, binary and 
ternary complex are described by macroscopic dissociation constants (K\, K 2 ). In addition, equations describing interdependencies of the different 
dissociation constants are given. K\ and K 2 as well as K J2 , Ki 3 , K 24 and K 3 4 are thermodynamic constants that describe the RNA-[Mg 2+ ] n complexes 
under equilibrium conditions. 



the ratio between dissociation constant (K D ) and RNA 
concentration and b is the maximum CSP at infinite 
Mg 2+ concentration. During the fitting procedure, the par- 
ameters a and b were allowed to adjust freely. 

CSP curves that showed biphasic behavior and could 
not be fitted by Equation (2) were fitted by a model 
assuming two different binding sites (Figure 1). Such a 
model is characterized by two macroscopic (K x , K 2 ) and 
four microscopic (K\, 2 , K\^, K 2 ,4, A^, 4 ) dissociation con- 
stants that connect four different RNA states (I, II, III, 
IV). Considering a particular imino group, two different 
states with only one Mg~ + ion bound (states II and III 
in Figure 1) can have different or equal (degenerated) 
chemical shifts. Such model predicts that biphasic 
curves only depend on the macroscopic dissociation 
constants K x and K 2 . Monophasic curves, in contrast, 
detect either macroscopic or microscopic dissociation 
constants, depending on which chemical shifts of the 
four states are degenerated for a particular imino group. 
Microscopic and macroscopic dissociation constants 
are related to each other by the equations depicted in 
Figure 1. 

Biphasic CSP curves can be fitted by Equations (3-5), 
if the two binary complexes have different chemical shifts. 



CSP = a*6 {BA 2)+b*e iAB) 
K X K 2 x [A] 



7 (AB) 



(BA2) 



[A] 2 

' (X^+^tAHA] 2 ) 



(3) 
(4) 

(5) 



Biphasic CSP curves can be fitted by Equation (5-7), if the 
complexes [AB] and [BA] are degenerated. 



CSP = a*6> ( BA2)+(6/2)* #(AB+BA) 



(AB+BA) 



K 2 x [A] 



{K,K 2 +K 2 [A]+[Af) 



(6) 



(7) 



Whether chemical shifts are degenerated in the binary 
complexes or not, influences only the parameter b in the 
Equations (3) and (6). In Equations (3) and (4), the con- 
stants b and Ki 3 cannot be fitted independently and can 
therefore be summarized to one constant. 

In Equations (3)-(7), [A] is the concentration of 



Mg -ions, 0(ab; an d #(ba) are the populations of the 
two different binary (RNA*Mg 2+ )-complexes and 6 f \B2j 
is the population of the ternary (Mg 2+ *RNA*Mg 2+ )- 
complex. Parameters a and b are fitting parameters that 
modulate the CSP amplitude. K\, K 2 are the macroscopic 
and Ki >2 , K 2 ^ 4 , K 3 4 the microscopic dissociation con- 
stants of the Mg 2+ binding model which assumes two 
Mg 2+ binding sites (Figure 1). During the fitting proced- 
ure, the parameters a, b, K x and K 2 were allowed to adjust 
freely. 

In some cases, three out of four RNA states have 
degenerated chemical shifts and only one state has a dif- 
fering chemical shift. For some of these cases, this results 
in monophasic CSP curves. For example, if state I has a 
differing chemical shift, the macroscopic dissociation 
constant K\ is detected. If state IV has a differing 
chemical shift from states I, II and III, the macroscopic 
dissociation constant K 2 is detected. In the special case of 
equal microscopic dissociation constants (A^ 1?2 , ^1,3, K 2A , 
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AT 3i4 ) the macroscopic dissociation constant K 2 is four 
times larger than the macroscopic dissociation constant 
K x (K 2 = A*K X ). 



MnCl 2 PRE experiments 

Paramagnetic relaxation enhancement (PRE) effects 
caused by paramagnetic Mn 2+ ions were analyzed for 
the observable imino groups within the 4U-hp2-wt 
RNA. ['H-^NJ-HSQC spectra of the imino region of a 
4U-hp2-wt RNA sample were recorded in NMR buffer 
in the presence of 5 mM MgCl 2 at MnCl 2 concentrations 
of 0 and 4 uM. Experiments were performed at a tempera- 
ture of 10°C and a B 0 field strength of 900 MHz (21.14T). 
The recycling delay was set to 3 s. Remaining acquisition 
parameters were identical to the MgCl 2 titration experi- 
ments. PRE effects were quantified by intensity and line 
width analysis. Line width analysis was performed for the 
indirect 15 N-dimension. Intensities and line widths were 
analyzed with Sparky 3.1.1.3. For the PRE analysis we 
plotted In (I[0 uM]/I[4 uM]) for each observable imino 
signal and compared the results with the change in line 
width (u[4 uM] - u[0 uM]). 



Co(NH 3 ) 3+ 6 NOESY experiment 

To identify and localize Mg 2+ outer sphere binding sites, a 
NOESY spectrum of a 4U-hp2-wt RNA sample in NMR 
buffer was recorded in the presence of 5 mM Co(NH 3 ) 3+ 6 . 
Measurements were performed at a temperature of 10°C 
and a B Q field strength of 600 MHz (14.09 T) with a mixing 
time T m of 150 ms. The spectrum was recorded with a 
spectral width of 15ppm in the indirect dimension and 
22 ppm in the direct dimension. A total of 896 increments 
were recorded in the /[ dimension and 2048 points in the t 2 - 
dimension. Each increment was recorded with 128 scans. 



RNA preparation 

Unlabeled and 15 N-labeled RNA oligonucleotides were 
synthesized by in vitro transcription with T7 RNA poly- 
merase from linearized plasmid DNA templates and 
purified as described previously (45). 15 N-labeled nucleo- 
tides were purchased from Silantes (Munich, Germany), 
unlabeled nucleotides were purchased from Sigma Aldrich 
(St. Louis, MO, USA). 

RNAs were desalted using Vivaspin 20 concentrators 
with a molecular weight cutoff (MWCO) of 3000 Da 
from Sartorius stedim biotech (Aubagne Cedex, France). 
RNA refolding was achieved by denaturing for 10 min at a 
concentration of 0.25 mM and a temperature of 95°C 
followed by 20-fold dilution with ice cold water and sub- 
sequent incubation at 0°C for another 30 min. RNAs were 
concentrated and exchanged to NMR buffer using 
Vivaspin 20 concentrators (MWCO = 3000 Da). The 
NMR buffer contained 15 mM K x H y P0 4 (pH 6.5), 
25 mM KC1, 90% H 2 0 and 10% D 2 0. Correct RNA 
folding was verified by native polyacrylamide gel 
electrophoresis. 



Native gel electrophoresis 

The 5x loading buffer was composed of 87% glycerole, 
0.1% (w/v) xylene cyanole FF and 0.1% (w/v) bromophenol 
blue. The 1 x running buffer contained 50 mM Tris-acetate 
and 50 mM sodium phosphate (pH 8.3). Gels were cast using 
15% acrylamide with an acrylamide:bisacrylamide ratio of 
37.5:1. Running conditions were set as follows: U = 70 V, 
P<0.5W for 4h, water cooling to prevent heating. Gels 
were stained with ethidium bromide. 

CD spectroscopy 

CD melting and refolding curves were recorded with a 
JASCO spectropolarimeter J-810 at a wavelength of 
258 nm. The RNA concentration was adjusted to 25 uM. 
Buffer conditions: 15mM K x H y (P0 4 ), 25 mM KC1, pH 
6.5. Melting curves were recorded with a temperature 
slope of l°C/min (4U-hp2-wt: 5-90°C; 4U-hp2-G14A- 
C25U-mutant: 2-80°C; 4U-hp2-A8C-mutant: 10-90°C) 
and refolding curves with a temperature slope of — 1°C/ 
min (4U-hp2-wt: 90-5°C; 4U-hp2-G14A-C25U-mutant: 
80-2°C; 4U-hp2-A8C-mutant: 90-10°C). CD melting and 
refolding curves were normalized according to Equation 
(9) and thermodynamic parameters (AH un f, AS un {, AG un f 
and T m ) were derived from Equations (10) to (12). 

Assuming the RNA to be either folded or unfolded 
(two-state model), the fraction of unfolded RNA can be 
expressed by 

[RNA fold ed(r)] ^ 
U(T)=X - [RNA total ] (8) 

The fraction of unfolded RNA a is varying with tem- 
perature. This temperature dependence of a can be 
determined from the temperature dependence of the ellip- 
ticity at a wavelength of 258 nm according to Equation (9) 

0(T) - £ U nfolded(7) 



a(T)=l 



&{o\ded(T) — ^unfolded (X) 



(9) 



In Equation (9), 6(T) is the measured temperature- 
dependent ellipticity, 0 un f o \d s d(T) and 0f o ided (T) are the 
ellipticities of the unfolded and the folded state of the 
RNA, respectively. 9 mfMed (T) and 9 fMed (T) are also 
slightly temperature-dependent and are therefore approxi- 
mated linearly. From the resulting a(T) dependence, the 
molar unfolding enthalpy AH un{ can be calculated (46) 
according to Equation (10). The melting point T m is 
defined as the temperature at which 50% of the RNA is 
unfolded (a = 0.5). In case of perfect sigmoidality, the 
slope of the a(T) curve is maximal at this temperature 



4RTi 



3a 
df 



(10) 



T=T m 



In Equation (11), R is the gas constant, T m the melting 
temperature, a the fraction of unfolded RNA and AH w ,f 
the molar unfolding enthalpy of the RNA. From AH un{ , 
the molar unfolding entropy A5 un f can be calculated ac- 
cording to Equation (11) 

AS unr = ^ (11) 
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From the Gibbs-Helmholtz Equation (12), the molar free 
energy AG M „/can be calculated, if one assumes A//„„yand 
AS un f to be temperature-independent 

AG unf = AH mf - TASurf (12) 

Mg 2+ -dependence of the RNA melting point T m 

The RNA melting point T m can be determined by CD 
spectroscopy. Assuming 7" m ([Mg]) to be linearly depend- 
ent on the fraction of RNA bound to Mg 2+ [Equation 
(13)], the free energy AG C d of Mg 2+ binding to the 
RNA can be calculated according to Equation (14): 

7m([M g 2+ ]) =r; M HAr m = ir Mg 

+A7 T m , max [RNA bound ]/[RNA t otai] 



r m ([M g 2+ ]) =r^+ 



at„ 



( [Mg 2+ ] 



e X p(AGc D /(RT))+[Mg 2+ ] 



jiloMg_ 



AT„ 



s [Mg 2+ ] 



(14) 



^CD+[Mg 2+ ] 



In Equation (14), T m is the melting point of the RNA, 
jnoMg j s me iting point in the absence of Mg 2+ ions, 
Ar m ,max is the maximal change in T m caused by Mg 2+ , 
[Mg 2+ ] is the Mg 2+ concentration, K CD is the RNA:Mg 2+ 
dissociation constant and AGcd is the free energy of Mg 2+ 
binding. 



bgaB reporter gene assay 

Three different translational fusion constructs have been 
used to study the temperature dependence of gene expres- 
sion in E. coli using the 6-galactosidase gene bgaB as 
reporter. Constructs were designed by introducing the 
temperature sensitive hairpin 2 of the 5'-UTR of the 
Salmonella agsA gene in front of the heat stable 
B-galactosidase gene bgaB utilizing an EcoRI restriction 
site (underlined in the sequences). 

The 4U-hp2-wt construct contains the second hairpin of 
the fourU (4U) RNA thermometer (41) in an optimized 
vector. Sequence: aGCGTTGAACTTTTGAATAGTGA 
TTCAGGAGGTTAATGATGGCAG_AATTC~/>g<zB. 
The 4U-hp2-A8C-mutant construct differs from the 
wild-type construct by an A to C mutation (highlighted 
boldly). Sequence: aGCGTTGACCTTTTGAATAGTGA 
TTCAGGAGGTTAATGATG GCA GAATTC ~6gfl£. 
In the 4U-hp2-G14A-C25U-mutant the G14-C25 base 
pair is exchanged to an A14-U25 base pair (highlighted 
boldly). Sequence: aGCGTTGAACTTTTAAATAGTGA 
TTTAGGAGGTTAATGATGGCAGAATTC-^afi. 
An amount of 25 ml LB-Amp (150ug/ml), pre-warmed in 
a shaker to 30° C were inoculated with 1.5 ml over night 
culture (5 ml grown at 30° C) of DH5a cells containing the 
plasmids. Cells were grown to an optical density (OD 60 o) 
of 0.5. After induction of transcripts with 25 ul L-arabin- 
ose (10%), 8ml of the culture were shifted to two 
pre-warmed, shaking flasks at different temperatures, re- 
spectively. After 30 min cells were used for the enzymatic 
assay as described before (47). To compare the different 



measurements the expression value for the wt at 30°C was 
set to 1, respectively. 



RESULTS AND DISCUSSION 

Hairpin stability of the fourU thermometer inversely 
correlates with gene expression 

In our previous study, we determined base-pair stabilities 
of hairpin 2 (4U-hp2-wt, Figure 2B), which represents the 
temperature responsive hairpin of the Salmonella fourU 
thermometer, in comparison with a stabilized point 
mutant (4U-hp2-A8C-mutant), in which the AG 
mismatch of the wt structure was replaced by a CG pair 
(Figure 2A). From this study, we concluded that the G14- 
C25 base pair is by far the most stable base pair in the 
thermometer hairpin and has a stabilizing effect on the 
secondary structure of the entire RNA helix (39). Thus, 
substitution of this stable GC base pair into a less stable 
AU base pair (Figure 2C) is predicted to result in a sig- 
nificant reduction of the thermal stability of the 4U-hp2 
RNA. CD unfolding curves recorded on the G14A-C25U 
double mutant indeed confirmed the expected behavior 
(Figure 3A). The GC to AU exchange resulted in a 
melting point T m of 35.9°C, which is 5.2°C lower than 
the T m value of the wild-type RNA and 16.3°C lower 
than the 4U-hp2-A8C mutant RNA (Figure 3A and B). 
T m values for the three RNAs under investigation are 
given in Table 1. 

Reversible melting and refolding of an RNA thermom- 
eter has recently been shown to be physiologically import- 
ant for fitness of a cyanobacterium (48). CD unfolding 
and refolding curves demonstrated a similar full reversibil- 
ity of the unfolding transition of all three fourU hairpins 
investigated in the current study since start and end 
points of unfolding and refolding curves are identical 
(Supplementary Figures S1-S3). However, unfolding and 
refolding curves slightly differ due to disturbance of the 
equilibrium caused by the continuous temperature change 
during the temperature scan. For the subsequent analysis, 
melting and refolding curves were averaged (the averaged 
curve will be referred to as unfolding curve from now on) 
to derive the unfolding curve under equilibrium condi- 
tions. Unfolding curves were normalized according to 
Equation (9) and a(T) values were determined assuming 
a two-state model (folded; unfolded) for RNA unfolding. 
The thermodynamic parameters AH un{ , AS un{ and AG unf 
were derived according to Equations (10-12) from the 
a(T) curves shown in Figure 3A. AH un{ , AS un f, AG un f 
and T m values for the three RNAs under investigation 
are given in Table 1. 

From the AH unt - and A5 un f values, the corresponding 
AG unt {T) curves can be calculated according to the 
Gibbs-Helmholtz equation. These AG unt {T) curves are 
illustrated for the wild-type RNA and the two mutants 
in Figure 3B. Remarkably, all three AG un f(T) curves 
exhibit an intersection point at T c = 81.8°C. At this tem- 
perature, the three RNAs possess the same stability 
AG unt {81.8°C) = -31.94kJ/mol. At T c , the unfolded 
RNA is almost exclusively populated since it is then 
much more stable than the folded conformation. The 
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(a) 



A8C mutant 



(c) 



wildtype 



G14A-C25U mutant 



4U motif 



A8C mutant 




Shine Dalgarno 



Figure 2. Sequence and secondary structure of hairpin 2 of the Salmonella fourU RNA thermometer and mutants used in this study. (A) Secondary 
structure of the 4U-hp2-A8C mutant RNA. (B) Secondary structure of the 4U-hp2-wt RNA. (C) Secondary structure of the 4U-hp2-G14A-C25U 
mutant RNA. 



existence of an intersection point of the AG unf {T) curves is 
a consequence of the linear enthalpy-entropy correlation 
of the global unfolding transition which is depicted in 
Figure 3C. The slope of the AH unt {AS un {) correlation is 
equal to m = 354.9 K and the ordinate intercept was 
determined to be yO = —31.95 ± 0.95kJ/mol. Linear 
enthalpy-entropy correlations for global unfolding transi- 
tions have been reported earlier by Searle and Williams 
(49) and by Petruska (50) for DNA duplexes. There is also 
a strong linear correlation between AG unf and AH unt - 
(Figure 3D) in the fourU RNA thermometer. The slope 
of this correlation is temperature-dependent and is 
illustrated for a temperature of 20°C. 

Insertion of the enthalpy-entropy correlation into the 
Gibbs-Helmholtz equation results in the following 
equation: 

AG unf = (m - T) x AS unf +jO (15) 

At the compensation temperature (T c = m), enthalpy and 
entropy contributions of the RNA-solvent interaction 
cancel each other and only the RNA-RNA interactions 
remain and equal yO. 

Helix unfolding can be approximated according to 
Searle and Williams (49) by the following equation 

AGunf ^ Ai/ S tacking — T X ASVotation (16) 

At T c , this leads to the following relation 

AG un f(T c ) = yO « A// Stacking — T x AS rotation (17) 



AG unt {T c ) is affected by the increase in entropy upon un- 
folding, mainly caused by the gain of rotational freedom 
in the RNA backbone ASVotation and the loss of stacking 
interactions between neighboring nucleobases A// Stac king- 
At the compensation temperature T c , the absolute 
value of the entropy term T*AS exceeds by far the 
absolute value of the enthalpy term A// Stacking . Thus, the 
RNA is almost completely unfolded at T c 
(AG unf (r c ) = yO = -31.95kJ/mol ± 0.95kJ/mol). 

In order to test whether the in vitro observations correl- 
ate with translational control in vivo, the 4U-hp2-wt RNA, 
the 4U-A8C RNA and the 4U-G14A-C25U RNA were 
analyzed in the bgaB reporter gene assay as described pre- 
viously (41). The 5'-UTRs were cloned 5'-upstream of the 
bgaB reporter gene downstream of an arabinose-inducible 
promoter and transformed into E. coli DH5a cells. 
Transcription was initiated with 0.01% (w/v) of L-arabin- 
ose and translation at various temperatures was followed 
by determination of P-galactosidase activity. Expression 
levels were normalized to the expression of the 
4U-hp2-wt sequence at 30°C. The gene expression curves 
were consistent with the in vitro data in that the least 
stable G14A-C25U mutant led to derepressed, perman- 
ently elevated reporter gene expression even at low tem- 
peratures (Figure 4). The wild-type 5'-UTR allowed 
intermediate gene expression at elevated temperatures 
whereas the stable A8C mutant blocked gene expression 
completely. The combined in vivo and in vitro results show 
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Figure 3. (A) Fraction of unfolded RNA ot calculated from the CD melting curves according to Equation (9) of the 4U-hp2-wt RNA (black line), 
4U-hp2-A8C mutant RNA (blue line) and the 4U-hp2-G14A-C25U double mutant RNA (red line) (B) AG un r(7") curves for the global unfolding of 
the 4U-hp2-wt RNA (black), A8C mutant RNA (blue) and the 4U-hp2-G14A-C25U double mutant RNA (red). AG unf <r) values were calculated 
from the corresponding A/Z^f and A5 un r values (Table 1) according to the Gibbs-Helmholtz equation. T m indicates the three different melting 
points (a = 0.5) of the three RNAs, respectively. T c = 354.9 K represents the compensation temperature at which all three RNAs have the same 
stability. (C) Enthalpy-entropy correlation for the global RNA unfolding determined from the three 4U-hp2 RNAs under investigation. (D) Free 
energy-entropy correlation of the global RNA unfolding (T = 20°C) determined for the three 4U-hp2 RNAs analyzed. Linear correlations in (C) and 
(D) were fitted according to the linear equation f = yO + mx. Fitting results are represented within the figures (C) and (D), respectively. 



Table 1. Thermodynamic parameters A// unf , A5 un f, AG mlf (T= 20°C) 
and T m for the global unfolding of the 4U-hp2-wt, 4U-hp2-A8C- 
mutant and 4U-hp2-G14A-C25U-mutant RNAs 





A-Z/unf 


ASunf 


AG un f 


T m 




(kJ/mol) 


[J/(mol*K)] 


(kJ/mol) 


(°C) 








(T= 20°C) 




Wild-type 


248.9 


792.0 


16.75 


41.1 


G14A-C25U mutant 


214.2 


692.9 


11.05 


35.9 


A8C mutant 


351.1 


1079 


34.81 


52.2 



that seemingly minor mutational alterations can have a 
very significant effect on RNA stability and functionality 
of a regulatory RNA element. By switching from the 
stable A8C sequence to the unstable G14A-C25U 
sequence, gene expression levels can be increased by 
factors up to 150. In addition, the wild-type RNA as 
well as the G14A-C25U RNA are able to induce gene 
expression in response to a temperature shift from 30° C 
to 42°C by factors of 8.6 (wt) and 3.9 (mut) (Figure 4). 
Induction factors of at least 3.5 are described for function- 
al thermoregulators (48,51-53) and an overall enhanced 
reporter gene expression for mutants with slightly 
destabilized structures has also been observed for the 
ibpA thermometer from E. coli (51). Therefore both 
ag^-5'-UTR-variants, wild-type and G14A-C25U, 



consisting only of the second hairpin of the natural 
occurring agsA thermometer, work well as control 
elements in the physiological range of living cells. We 
observed correlation between hairpin stability determined 
in vitro and gene repression levels measured in vivo. 
Interestingly, gene expression levels start to increase at 
temperatures significantly lower than the melting point 
determined in vitro. Apparently, translation initiation 
already occurs even at temperatures where only a minor 
fraction of RNA thermometers are unfolded. 

Curvatures of the temperature dependent gene expres- 
sion curves in vivo (Figure 4) differ from the CD unfolding 
curves measured in vitro (Figure 3A and B). Additional 
effects to the accessibility of the SD sequence seem to 
modify the gene expression in vivo. Other possible 
factors like the temperature dependence of transcription, 
RNA degradation or translation initiation may have add- 
itional modulating effects that superimpose the initial 
RNA thermometer effect. 

Mg 2+ -dependence of the melting point 

One well-known factor that contributes to folding of 
complex RNA structures is Mg 2+ . In addition, based on 
the sequence of the fourU thermometer, a Mg 2+ binding 
site can be predicted (32). To examine whether this ion 
plays a structural and functional role in the fourU therm- 
ometer, we monitored the Mg 2+ dependence of the melting 



Nucleic Acids Research, 2011, Vol. 39, No. 18 8265 



protein expression levels 





▲ G14A-C25U 

♦ WT 

A A8C 

* * 


A 

A A 
A * 

A 

* 

♦ ♦ ♦ ♦ 





26 28 30 32 34 36 38 40 

temperature [°C] 



Figure 4. Temperature dependent expression of the agsA-(4U-hp2)- 
bgaB fusion in DH5a E. coli cells. Following RNA constructs were 
investigated: 4U-hp2-wt (filled diamonds), 4U-hp2-A8C-mutant (open 
triangles), 4U-hp2-G14A-C25U (filled triangles). Expression levels were 
referenced to the respective expression level of the 4U-hp2-wt at 30°C 
which was set to 1. 
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Figure 5. T m melting point dependence of the 4U-hp2-wt RNA on the 
MgCl 2 concentration. Data points are given as filled diamonds while 
the fit according to Equation (14) is illustrated as continuous line. 



Table 2. Fitting results of the 7' m (Mg 2+ ) dependence of the 
4U-hp2-wt RNA according to Equation (14) 



4U-hp2-wt 7^ oM s 


A J m, max 


^CD 


AG CD 


(°C) 


(°Q 


(uM) 


(kj/mol) 


41.77 ± 0.22 


16.93 ± 0.25 


899 ;l ± 47.9 


-18.860 ± 0.14 



^Kcd was calculated for a temperature of 50.24°C. Errors are derived 
from the fit. 



point T m of the 4U-hp2-wt RNA by CD spectroscopy. 
The T m (Mg 2+ ) dependence was fitted by Equation (14) 
(Figure 5) and AG CD , Ar mmax and T"° Mg values were 
derived (Table 2). 

The 7" m (Mg 2+ ) dependence can be described by a model 
which assumes only a single binding site. NMR experi- 
ments, however, provide evidence for two Mg 2+ binding 
sites (see below). Since Mg 2+ binding to these sites is co- 
operative (see below), the 4U-hp2-wt RNA exists most 
prominently in either the free form or the ternary 
complex ([Mg 2+ *RNA*Mg 2+ ]) but rarely in one of the 
binary complexes ([Mg 2+ *RNA] or [RNA*Mg 2+ ]). Thus, 
the r m (Mg + ) curve exhibits hyperbolic shape and the 
AG C d value determined from the r m (Mg 2+ )-dependence 



represents the free binding energy between free RNA and 
ternary [Mg 2+ *RNA*Mg + ] complex. In addition, we can 
show that the temperature dependence of T m is caused by 
specific and not by diffuse Mg 2+ binding (Supplementary 
Data: Evidence for defined Mg 2+ binding). The fact that 
the temperature dependence is dominated by defined 
Mg 2+ binding does not rule out the possibility of addition- 
al diffuse binding. However, according to our data the 
effects of diffuse Mg 2+ binding are small in comparison 
to the effects caused by specific binding, at least in the 
concentration range observed here (0-7 mM). 

The Mg 2+ dependence of T m observed in vitro indicates 
a potential effect of fluctuating Mg 2+ concentrations 
within the cell on expression levels of genes under the 
control of the fourU RNA thermometer. According to 
the literature, the free intracellular Mg 2+ concentration 
may vary between 1 and 2mM (14,18). Such a shift in 
Mg 2+ concentration causes a shift of the T m value of the 
4U-hp2-wt RNA by 2.8°C (Figure 5). Thus, varying Mg 2+ 
concentration will most probably have a modulating 
temperature-independent effect on RNA-thermometer- 
controlled gene expression in vivo. 

Localization of the Mg 2+ binding sites 

We applied different NMR methods to identify Mg 2+ 
binding sites in the 4U-hp2-wt RNA. First, a NOESY 
experiment of a 4U-hp2-wt RNA sample in the presence 
of [Co(NH 3 ) 6 ] 3+ ions was performed to identify and 
localize [Co(NH 3 ) 6 ] 3+ binding sites. [Co(NH 3 ) 6 ] 3+ ions 
serve as mimic for hexahydrated Mg 2+ ions and therefore 
map binding sites for outer sphere coordinated Mg 2+ ions 
(28,29,32). Second, we analyzed the PRE effects of the 
imino signals occurring upon the addition of paramagnet- 
ic Mn 2+ ions to the 4U-hp2-wt RNA sample (28,29). The 
fact that substoichiometric amounts of Mn 2+ are sufficient 
to cause significant line broadening of imino signals in 
close proximity to Mg 2+ binding sites reveals that Mg 2+ 
ions are in rapid exchange with free divalent ions. The use 
of Mn 2+ as probe for Mg 2+ is supported by the observa- 
tion that 76% of the intracellular Mg 2+ ions in prokary- 
otes (E. coli) can be substituted by Mn 2+ ions without 
severely impairing viability of the cells (18). Third, we 
followed CSPs of the imino signals during Mg 2+ titration 
experiments (28,29). Fitting of the obtained CSP curves 
was performed to derive K D values for Mg 2+ binding. By 
taking the results of these three complementary techniques 
into account, it became possible to unravel the different 
effects on how Mg 2+ ions interact with the 4U-hp2-wt 
RNA. Combination of the information yields a consistent 
overall picture on the effects of Mg 2+ ions on the 
4U-hp2-wt RNA. 

The imino region of the 2D- 1 H- 1 H NOESY spectrum of 
the 4U-hp2-wt RNA in the presence of 5mM 
[Co(NH 3 ) 6 ] 3+ ions is given in Figure 6. Protons from 
[Co(NH 3 ) 6 ] 3+ exhibit strong cross peaks to the imino 
protons of the nucleobases Ull, U12, G27 and/or G28 
and weak cross peaks to the imino protons of the 
adjacent nucleobases U10 and U13. Thus, a binding site 
for [Co(NH 3 ) 6 ] 3+ is located in close proximity to the 
nucleobases Ull, U12, G27 and G28, most probably in 
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the major groove of the RNA helix very similar to the 
binding site described for the structure of the P5b 
hairpin of the group I intron by Kieft and Tinoco (32). 

PRE effects of the imino signals caused by the addition 
of Mn 2+ ions were analyzed in order to reveal divalent 
cation binding sites. For that purpose, [ 1 H, 15 N|-HSQCs 
of the imino region of the 4U-hp2-wt RNA sample were 
recorded and line shapes and intensities of the imino 




15 14 13 12 11 10 



OJ 1 - 1 H (ppm) 

Figure 6. Section from the 2D-'H NOESY experiment performed on 
a 0.9-mM 4U-hp2-wt RNA sample in the presence of 5 mM 
[(Co(NH 3 ) 6 ] 3+ . Assignments of cross signals between protons of the 
[Co(NH 3 ) 6 ] and the imino protons of the RNA are given within 
the spectrum. On top of the NOESY spectrum, the corresponding 
imino region of the 'H-1D NMR spectrum is given. Buffer conditions: 
15mM K x H y P0 4 (pH 6.5), 25 mM KC1, 90% H z O and 10% D 2 0. 



signals were analyzed to quantify the PRE effects. 
Changes in line widths and signal intensities caused by 
the addition of 4uM Mn 2+ are given in Figure 7. Large 
PRE effects were observed for the nucleobases Ull and 
U12 within the RNA helix. In addition, the nucleobases 
U13, G14 and U33 also show significant PRE effects, 
while the imino signal of nucleobase U10 exhibits only a 
very weak PRE effect. Interestingly the G27/G28 double 
signal only shows a moderate PRE effect in comparison to 
Ull, U12 and U13 which is consistent with binding of 
[Co(NH 3 ) 6 ] 3+ in the major groove. 

Taken together, the PRE analysis reveals the existence 
of a Mn 2+ binding site within the 4U-hp2-wt RNA which 
is located in close proximity to the nucleobases Ull, U12 
and U13. This Mn + binding site represents a strong indi- 
cation for the existence of a Mg 2+ binding site at exactly 
the same position within the 4U-hp2-wt RNA. 
Furthermore, PRE analysis and [Co(NH 3 ) 6 ] 3+ NOESY 
data reveal binding sites in the same region of the RNA 
leading to the conclusion that there is one specific outer 
sphere Mg 2+ binding site within the fourU motif of the 
4U-hp2-wt RNA. Moderate PRE effects of the 
nucleobases G27/G28 are in agreement with localization 
of this binding site within the major groove. 

Finally, we titrated 4U-hp2-wt RNA with Mg 2+ ions 
and analyzed the occurring CSPs of the imino group 
region within the ['lI, 15 N]-HSQC spectra. CSP([Mg 2+ ]) 
curves were fitted according to Equation (2) and K D 
values were derived. Figure 8 shows Mg 2+ -dependent 
CSPs of the imino groups with their respective fits 
(continuous lines). In Figure 8D, the shifting signals of 
the imino groups U5 and U13 within the ['H- 15 N]- 
HSQC during the titration of the RNA with MgCl 2 
are shown exemplarily. Fitting results are summarized 
in Supplementary Table SI. Interestingly, different 



MnCI 2 PRE of the wildtype RNA 




G28 

nucleotide 

Figure 7. PRE of the imino resonances of the 4U-hp2-wt RNA observed upon the addition of 4 uM MnCl 2 . Changes in line widths Av [Hz] in the 
15 N-dimension are illustrated as filled triangles (continuous line) with the corresponding scale on the left hand side of the diagram and changes in 
intensity (ln(I[0 uM]/I[4 uM]) [a.u.]) are illustrated as open squares (dotted line) with the corresponding scale on the right hand side. U32 and G31 do 
not give rise to any observable imino signals and the imino signals of G30 and U4 are too weak to be analyzed. 
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Figure 8. (A-C) CSPs of the imino resonances of the nucleobases U13, U23, G14, U24, U4, U5 and G6 plotted against the [Mg 2+ ]/[RNA] ratio and 
fitted according to Equation (2) (continuous lines). CSPs of the imino signal U12 (B) were fitted by Equation (3) (D) Mg 2+ -dependent changes (CSPs) 
of the imino groups U13 and U5 within the ['H, 15 N]-HSQC spectrum (titration steps: see 'Materials and Methods' section). 



K diss values are obtained for different imino signals. 
Lowest K diss values are found for U4 (0.85 ± 0.18) 
and G14 (1.17 ± 0.12 mM). Nucleotides U5, G6 and 
U13 exhibit K diss values of around 2.34 mM and the nu- 
cleotides U23 and U24 show K diss values of ~3.56mM. 
U13 shows by far the largest CSP max value 
(0.250 ± 0.008 ppm). 

Multiple conformations of a [RNA*Mg 2+ ]-complex 
assuming only a single Mg 2+ binding site fail to explain 
the observed diversity in ^ diss values as explained in the 
Supplementary Data (Mg 2+ binding model for one 
binding site but different conformations of the 
[RNA*Mg 2+ ] complex). 

In addition, the CSP curve of U12 shows a distinct 
biphasic behavior. We argue that this behavior can only 
be explained by two different Mg 2+ binding sites. The 
model describing the binding properties of an RNA 
molecule with two different Mg + binding sites is given 
in Figure 1. The dissociation constants derived from the 
biphasic fit (CSP curve of U12) can be assigned to the 
macroscopic dissociation constants K\ (3.58 mM) and 
K 2 (2.06 mM). Since the observed CSP curves cannot 
be explained by negative cooperativity (negative 
cooperativity would lead to negative microscopic dissoci- 
ation constants), cooperativity has to be positive. In such 
model, the smallest dissociation constants derived by 
monophasic fitting have to be purely monophasic and cor- 
respond to one of the microscopic dissociation constants 
K 2A and K 3A . Thus, K 2A was set to 1.17 ± 0.12mM (G14) 



and K 3A to 0.85 ± 0.18 mM (U4) as derived from 
monophasic fitting of the CSP curves (Figure 8 and 
Supplementary Table SI). However, an assignment of 
the microscopic dissociation constant K 2A (1.17mM) or 
the microscopic dissociation constant K 3A (0.89 mM) to 
the outer sphere Mg 2+ binding site within the fourU 
motif is not possible. According to the equations connect- 
ing the different dissociation constants (Figure 1), the 
microscopic dissociation constants K X 2 and K\j could 
be derived (Table 3). Free energy values for Mg 2+ 
binding AGki and AGk2 at a temperature of 10°C as 
well as the free energy for ternary [Mg 2+ *RNA*Mg 2+ ] 
complex formation AG ges are given in Table 3. 

The macroscopic dissociation constant K\ is notably 
larger than K 2 although formation of the binary 
complex is statistically favored over the formation of the 
ternary complex. The resulting microscopic dissociation 
constants connecting the free RNA and the binary 
complexes K X 2 (6.31 mM) and K X3 (8.26 mM) are much 
larger than the corresponding dissociation constants K 2A 
(1.17mM) and K 34 (0.89 mM) connecting the binary 
complexes with the ternary complex. In other words, 
binding of one Mg 2+ ion promotes binding of the 
second Mg 2+ ion (positive cooperativity). According to 
the determined dissociation constants (Table 3), we deter- 
mine a Hill coefficient b = 1.47 (Supplementary 
Figure S4). This intermediate positive cooperativity is in 
line with positive cooperativity values for Mg~ + binding to 
tRNA phe reported in the literature (b = 1.7) (54). 
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Table 3. Thermodynamic parameters of Mg 2+ binding as well as 
macroscopic and microscopic dissociation constants of the 
two-binding-site-model (Figure 1) 



K Jlss (mM) 




3.58 ± 0.67 mM 


K 2 


2.06 ± 0.66 mM 


K 2 ,4 


1.17 ± 0.12mM 


^3,4 


0.85 ± 0.12mM 


K la 


6.31 ± 2.43 mM 




8.68 ± 3.71 mM 


AG K1 (10°C) 


-13.262 ± 0.446 kJ/mol 


AG K2 (10°C) 


-14.556 ± 0.780 kJ/mol 


AG ges (10°C) 


-27.818 ± 0.899 kJ/mol 


AG CD (50.24°C) 


-18.860 ± 0.143 kJ/mol 


AG ges (37°C) 


-21.807 ± 0.443 kJ/mol 


A// ges 


-90.855 ± 6.469 kJ/mol 


7*A5 ges (10°C) 


-63.037 ± 6.404 kJ/mol 



AG K i, AG K 2 an d AG ges were calculated for a temperature of 10°C. 
AGcd (50.24°C) was derived from CD melting studies. AH ges , 
r*AS g es and AG ges (37°C) were calculated from the difference of 
AGges (10°C) and AG CD (50.24°C). Errors in K u K 2 , K 2A , K 3A and 
^cd (Table 2) are statistical errors resulting from the fit. Errors in 
K U2 , K h3 , AG K1 (10°C), AG K2 (10°C), AG ges (10°C), AG CD 
(50.24°C), A//ges and r*zlS g es (10°C) are Gaussian errors calculated 
by error propagation. 

A description how two Mg 2+ binding sites may lead to 
the observed variety of monophasic CSP curves is given 
in the Supplementary Data (Qualitative explanation of the 
observed CSP curves using a two binding site model). 

From the Mg 2+ dependence of the melting point T m 
determined by CD spectroscopy (Figure 5 and Table 2), 
the free energy for Mg 2+ binding AG CD at a temperature 
of 50.24°C was determined. Assuming that Mg 2+ only 
binds to the folded but not to the unfolded RNA A6 CD 
is equal to AG ges . From the difference of AG ges 
(T=10°C) and AG CD (T = 50.24° C), AH ges 
and T*AS ges (10°C) were calculated according to the 
Gibbs-Helmholtz equation assuming that AH ges and 
ASges are not temperature-dependent. Results are given 
in Table 3. 

The T*AS„ es value (-63.037 ± 6.404 kJ/mol) indicates 
that the Mg + -bound state exhibits a higher degree of 
order than the free RNA which appears reasonable since 
in the ternary [Mg 2+ *RNA*Mg ] complex, rotational 
and translational degrees of freedom of the hexahydrated 
Mg 2+ ions are reduced. In particular, binding of two 
hexahydrated Mg 2+ ions to the RNA molecule is 
accompanied by the loss of three rotational and three 
translational degrees of freedom per hexahydrated Mg 2+ 
ion while the loss of one degree of freedom reduces the 
entropy by Vi R. However, this accounts only for 50 J/ 
(mol*K) (2*6/2*R) which is far less than the observed 
change in entropy {AS ge!l = 223 ± 23 J/(mol*K)). Thus, 
additional entropy reducing effects have to be present. 
Such effects might be caused by an increased conform- 
ational restriction of the RNA itself or by a higher 
degree of order in the water molecules involved. This en- 
ergetically unfavorable entropy term is overcompensated 
by stabilizing AH ges contributions (-90.855 ± 6.469 kJ/ 
mol) which most probably result from energetically favor- 
able Coulomb interactions. 



Using AGcd (50.24°C) and AG ges (10°C) values 
(Table 3), the free energy for Mg 2+ binding AG ges (37°C) 
calculates to -21.807 ± 0.443 kJ/mol. Such value indi- 
cates that most probably not more than two Mg 2+ 
binding sites are involved since three or more Mg 2+ 
binding sites would lead to a lower AG ees value or to 
K diss values higher than those reported for diffuse 
binding (55). It may be noticed that the average AG ges 
value that we obtained per magnesium binding site 
(— 10.9 kJ /mol) is in surprisingly good agreement with 
the average value (~-llkJ/mol) that can be derived 
from the difference in concentration of bound and free 
magnesium in vivo (13,14). 

In summary, these results indicate the existence of two 
Mg 2+ binding sites within the 4U-hp2-wt RNA. At least 
one of these binding sites is outer sphere coordinated 
and is located in proximity to the nucleobases Ull, 
U12, G27 and G28 most probably within the major 
groove of the RNA. Although the existence of a second 
Mg 2+ binding site is certain, it could not be localized 
unambiguously. PRE experiments suggest that it is near 
the nucleobase U33 (Figure 7). The fact that only one 
binding site could be located with [Co(NH 3 ) 6 ] 3+ suggests 
that the second Mg 2+ binding site might be inner sphere 
coordinated. However, it is also possible that the two 
binding sites are in such close proximity that they 
cannot be identified as two independent binding sites. 

Using a combination of in vitro and in vivo experiments 
we show that fourU RNA thermometers are optimized in 
their sequence and Mg 2+ binding properties. Two Mg 2+ 
binding sites with positive cooperativity could be 
identified within the 4U-hp2-wt RNA, one of them as 
outer sphere binding site directly within the fourU motif 
as predicted from the secondary structure of this RNA 
(32). Mg 2+ concentration has a modulating effect on the 
melting point of the fourU RNA thermometer and there- 
fore most probably influences gene expression in vivo. 
From the in vivo experiments it became obvious that the 
onset of gene expression already occurs at temperatures 
significantly below the melting point. Apparently, a small 
fraction of RNA thermometers in the open conformation 
seems to be sufficient to efficiently permit translation ini- 
tiation. In addition, the in vivo experiments reveal that 
already small changes in the temperature dependent sta- 
bility (T m ) of the RNA thermometer lead to tremendous 
effects on gene expression of the downstream gene. 
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